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Abstract. The present lack of a stable method to compare persistent homology groups with torsion 
is a relevant problem in current research about Persistent Homology and its applications in Pattern 
Recognition. In this paper we introduce a pseudo-distance dx that represents a possible solution to 
this problem. Indeed, dx is a pseudo-distance between multidimensional persistent homology groups 
with coefficients in an Abelian group, hence possibly having torsion. Our main theorem proves the 
stability of the new pseudo-distance with respect to the change of the filtering function, expressed both 
with respect to the max-norm and to the natural pseudo-distance between topological spaces endowed 
with R"-vaIued filtering functions. Furthermore, we prove a result showing the relationship between 
dx and the matching distance in the 1-dimensional case, when the homology coefficients are taken in 
a field and hence the comparison can be made. 



Introduction 

The extension of the theory of Persistent Homology to the case of persistent homology groups with 
torsion is an important open problem. Indeed, the lack of a complete representation similar to the one 
given by persistent diagrams is a relevant obstacle to this development. However, this fact does not 
prevent us from using persistent homology groups with torsion for shape comparison. In order to show 
this, we introduce here a new pseudo-distance dx between persistent homology groups with coefficients 
in an Abelian group (hence possibly having torsion) , proving that dx is stable with respect to changes of 
the filtering functions, measured by the max-norm. We also express this stability as a lower bound for 
the natural pseudo-distance between topological spaces endowed with M"-valued filtering functions. This 
approach opens the way to the use of persistent homology groups in concrete applications, extending 
some ideas developed in [22 for size homotopy groups to persistent homology. 

Now, let us illustrate how persistent homology groups with torsion can appear in applications. As 
an example, let us imagine to have to study all the possible ways of grasping an object by a robotic 
hand endowed with m fingers. Obviously, not all grasps can be considered equivalent. For example, 
some of them can be difficult to realize in practice, due to the position of the object and the technical 
capability of the hand. Therefore, it is reasonable to associate each grasp with a "cost" , represented 
by a real number. In one of the simplest cases, we can consider two fingers grasping an object whose 
surface S is diffeomorphic to a spherical surface, in presence of friction. Each grasp can be seen as an 
unordered pair {p, q} of points on S, and we can consider the set G of all possible grasps (endowed with 
the Hausdorff distance between unordered pairs). If the cost of the grasp is given by a value tpg{{p, q}), 
we are naturally led to the problem of comparing the similarity of the grasps possible for two different 
objects with surfaces S and S' , with respect to the cost functions tpg : C/ — >• K, (fg> : C/' — > K. This 
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can be done by computing the (singular) persistent homology groups of each topological space Q with 
respect to the filtering function (pg. In other words, we can consider the persistent homology groups 
Hf''^^\u,v) given by the equivalence classes of cycles in Hk{Gv) that contain at least one cycle in 
where Qt represents the set of points of Q at which tpg takes a value not greater than t. It is interesting 
to note that, in general, this procedure leads to groups with torsion, when the homology coefficients are 
the integer numbers. 

To understand this, let us consider a particularly simple example. Let us assume that the cost of 
the grasp is given by the value 'Pgi{p,q}) = —\\p — q\\, so relating the cost of a grasp to the level of 
closure of the robotic hand. For the sake of simplicity, let us examine the case S = = {{x,y,z) G 
: a:^ + + = 1}. In this case, the group h[^''^^\u,v) has torsion, when u < v < and u is 
close enough to zero. Indeed, we can consider the continuous function / taking each unordered pair 
{Pi ?} G Gv to the direction of the line I through p and q (note that {p, q} S Gv implies p ^ q, so that I 
is well-defined). This direction can be seen as a point in the real projective plane MP^. Let us observe 
that, for V close enough to zero, the function f '■ Gv ^ is a homotopy equivalence (the homotopy 
inverse of / is the function g : MP^ Gv that takes each direction w to the unordered pair {p, g}, 
where p and q are the points at which the line through (0,0,0), having direction w, meets S"^). Hence 
/ induces an isomorphism between the singular homology groups of Gv and MP . As a consequence, 
since i/i(IRP^) has torsion, the persistent homology group Hf''^'^\u^v) has torsion, too, when u is close 
enough to zero. For an introduction to Topological Robotics we refer the interested reader to |21| . 

Multidimensional persistence. Persistent homology has turned out to be a key mathematical method 
for studying the topology of data, with applications in an increasing number of fields, ranging from shape 
description (e.g., [51 IHl [Ml US] ) to data simplification |20] and hole detection in sensor networks [Tl]. 
Recent surveys on the topic include [TH Uni [HI [2^ . Persistent homology describes topological events 
occurring through the filtration of a topological space X (e.g., creation, merging, cancellation of con- 
nected components, tunnels, voids). Filtrations are usually expressed by real functions (/s : X M 
called filtering Junctions. The main idea underlying this approach is that the most important piece 
of information enclosed in geometrical data is usually the one that is "persistenf" with respect to the 
defining parameters. 

Until recently, research on persistence has mainly focused on the use of scalar functions for describ- 
ing filtrations. The extent to which this theory can be generalized to a situation in which two or more 
functions characterize the data is currently under investigation [1] [51 [31 [H [7] . This generalization to 
vector-valued functions is usually known as the Multidimensional Persistence Theory, where the adjec- 
tive multidimensional refers to the fact that filtering functions are vector-valued, and has no connections 
with the dimensionality of the space under study. The use of vector-valued filtering functions in this 
context enables the analysis of richer data structures. 

An important topic in current research about multidimensional persistent homology is the stability 
problem. In plain words, we need to determine how the computation of invariants in this theory is 
affected by the unavoidable presence of noise and approximation errors. Indeed, it is clear that any 
data acquisition is subject to perturbations and, if persistent homology were not stable, then distinct 
computational investigations of the same object could produce completely different results. Obviously, 
this would make it impossible to use such a mathematical theory in real applications. 

Up to our knowledge, no theoretical result is available for the metric comparison of persistent homol- 
ogy groups in presence of torsion, at the time we are writing. 

Prior works. The problem of stability in persistent homology has been studied by Cohen-Steiner, 
Edelsbrunner and Harer in [llj for scalar filtering functions. By using a descriptor called a persistence 
diagram, they prove that persistent Betti numbers are stable under perturbations of filtering functions 
with respect to the max-norm, provided that the considered filtering functions are tame. The same 
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problem is studied in [T2] for tame Lipschitz functions. In [10], Chazal et al. use the concept of 
persistence module and prove stability under the assumption that it is finite-dimensional. The problem 
of stability for scalar filtering functions is also approached in |13) , where it is solved by assuming that 
the considered filtering functions are no more than continuous, but only for the 0th homology. 

Multidimensional persistence was firstly investigated in [22^ as regards homotopy groups, and by 
Carlsson and Zomorodian in 6 as regards homology modules. In this context, the first stability result 
has been obtained for the 0th homology in [1 : A distance between the 0th persistent Betti numbers, also 
called size functions, has been introduced and proven to be stable under perturbations of continuous 
vector- valued filtering functions. Such a result has been partially extended in [2 for all homology 
degrees, under the restrictive assumption that the vector-valued filtering functions are max-tame. 

Contributions. In this paper we introduce a new pseudo-distance dr between persistent homology 



groups, applicable also in the case that these groups have torsion (Theorem 2.2). Moreover, we prove 



that dT is stable with respect to perturbation of the filtering function, measured by the max-norm 



(Corollary 2.101. This result is a consequence of a new lower bound for the natural pseudo-distance 5 



(Theorem 2.8), which also implies an easily computable lower bound for 5 (Corollary 2.9). Furthermore, 
we investigate the link between dj- and the matching distance d„iatch between persistent diagrams, in 
the particular case that the filtering functions take values in M and the coefficients are taken in a field, 
so that also dmatch can be applied (Theorem |2.4[ ) . 

The structure of the paper is as follows. In Section [l] we recall the definition of persistent homology 
group and show some examples. In Section [2] the definition of dx is given and our results are proven. 

1. Some definitions and properties in multidimensional persistence 

In this paper, the following relations ^ and -< are defined in M": for u = (wi,...,u„) and v = 
(wi,...,w„), we say u < v (resp. u ^ v) \i and only if Ui < Vi (resp. Ui < Vi) for every index 
i ^ 1, . . . ,n. Moreover, R" is endowed with the usual max-norm: ||(wi,M2, . . . — maxi<i<„ \ui\. 

We shall use the following notations: A+ will be the open set {(m, v) € M" x M" : u ^ v}. For every 
n-tuple u = {ui, . . . , Un) £ K" and for every function : X ^ K", we shall denote by X{if ^ u) the set 
{x X : ipi{x) < Ui, i = 1, . . . , n}. If X is a compact topological space and if : X ~> M" is a continuous 
function, we shall set ||<p||oo — T^so^xex \\0{x)\\oo- 

Now we can recall the definition of multidimensional persistent homology group. 

Definition 1.1 (Persistent homology group). Let fc g Z. Let X be a topological space, and ip : X ^ M" 
a continuous function. Let i^ : Hk{X{ip ^ u)) — !■ Hk{X{ip ^ v)) be the homomorphism induced by the 
inclusion map i : X{lp < u) X{0 < v) with u < v, where Hk denotes the fcth homology group with 
coefficients in an Abelian group G. If u ^ v, the image of is called the multidimensional kth persistent 
homology group of {X, ip) at {u, v). We shall denote it by the symbol h\^''^\u, v). The function h\^''^'^ 
is called the multidimensional kth persistent homology group of {X, ip). 

In other words, the group H^j^'"^^ (u, v) contains all and only the homology classes of cycles born 
before or at u and still alive at v. We observe that, up to Definition [LT] each persistent homology group 
is a function h\^''^\ taking each pair [u, v) e A+ to a group. 

In the rest of this paper, the superscript arrow will be omitted when the filtering function takes values 
in M. 

Remark 1.2. The choice of the particular homology we use is not essential in this context, in the sense 
that the results proven in this paper do not depend on this choice. However, while we will not further 
discuss this point, we recall that some important properties of persistent homology groups depend on 
the homology that is considered. As an example, the use of Cech homology with real coefficients allows 
us to assume that the function : A+ — > N U {oo} that takes each pair [u, v) E A+ to the rank of 
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the multidimensional kth persistent homology group of {X, (p) at (u, v) , is right-continuous in both its 
variables (cf. 0). 

Example 1.3. Let MP'^ be the real projective plane, represented by the unordered pairs {p, — p} of 
opposite points of 5*^ = {{x,y,z) £ : + y"^ + = 1}. Let us consider the filtering functions 
tp,ilj : MP^ -> E, defined by setting ip{{{x,y, z),—{x,y, z)}) = \z\ and 'tp{{{x,y, z),—{x,y, z)}) — 2\z\. 
Let us also consider the filtering function x ■ S'^ ^ defined by setting z) — \z\. For each 

{u,v) G A^, the persistent homology groups in degree 1 computed at (it,f) of if, x ^^'^ 4' ^^e displayed 
in Figure [T] from left to right. The coefficients of the homology are taken in Z. 




Figure 1. Three examples of persistent homology groups with integer coefficients, in 
degree 1. In the first and second case the topological space is the projective plane, 
represented by the unordered pairs {p, — p} of opposite points of S'^ , while in the third 
case the topological space is 5^. The filtering functions are (/3({(x, y, z), — (a;, y, 2;)}) = 
|z|, ip{{{x,y, z),—{x,y, z)}) = 2\z\ and z) = |z|, respectively. In each region, 

the persistent homology group at the points of that region is displayed. The symbol 
denotes the trivial group. 

Since we do not take coefficients in a field, the homology groups we consider have torsion, in the 
general case. As a consequence, they cannot be described by persistence diagrams (cf. [2]). It follows 
that the classical matching distance cannot be applied, if we are interested in the information given by 
the torsion part of our groups. However, we can define a pseudo-distance that does not suffer from this 
limitation. We shall do this in the next section. 

2. The pseudo-distance cLt and the proof of its stability 

In the following. Ah will represent the collection of all Abelian groups, considered up to isomorphisms. 

Definition 2.1. Let X, Y be two topological spaces, and Cp : X ^ M", t\j -.Y ^ M" two continuous 

filtering functions. Let i/^^''^'' : A+ Ah and Hj^'^'^ : A+ Ah the multidimensional persistent k- 

homology groups associated with the pairs {X, 0) and (Y, i/j), respectively. We assume that the homology 
coefficients are taken in an Abelian group. Let us consider the set i? of all e > such that, setting 
e~ (e, . . . , e) G M", the following statements hold for each {u, v) £ A+: 

(1) a surjective homomorphism from a subgroup of Hi^'^'^ {u, v) onto Hj^^^^ {u-e,v + e)} exists; 

(2) a surjective homomorphism from a subgroup of (^^"^^ (u, v) onto h[^^^^ iu-e,v + e)} exists. 

We define dx iH\^''^\ ^j^'^^j equal to inf if i? is not empty, and equal to 00 otherwise. 
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Before proceeding, we recall that the term "pseudo-distance" means that the considered function 
verifies all the properties of a distance, with the possible exception of the one assuring that two elements 
having a vanishing distance must coincide. The term "extended" refers to the possibility that the 
function takes the value oo. 

Theorem 2.2. The function dx is an extended pseudo-distance. 

Proof. The following statements hold: 

i) : dr is non- negative by definition; 

ii) : The equality dx {^h\^''^\ h\^''^^^ = is obtained by taking e = and considering the identical 

homomorphism id : H^^f^''^^ (u, v) — > H^^f^''^^ (u ~ e^v + e); 
Hi): dx is symmetrical by definition; 

iv): Assume that dr (^hI^^'^\ Hj^'"^^^ = ei and dr (^hI^''^\ Hjf'^^^ = 63. Let tj be an arbitrarily 
small positive real number. The definition of dx implies that for each (u, v) € A"*" a surjective 
homomorphism / from a subgroup F of H^''^\u^v) onto Hj^'^\u — {ei + ff),v + {ci + ff)), 
and a surjective homomorphism g from a subgroup G of Hj^'^'^ (u — (ei + 77) , w + (ei + ffj) onto 
h\^'^\u — (e*! + jf) — {t2 + rfjjV + (e*! +ff) + (£2 + v)) exist. Now, we can consider the subgroup 
F ^ f-^{G) C h[^''^\u,v). Obviously, the function g\a o fp is a surjective homomorphism 
from the subgroup F of h\^''^\u, v) onto Hjf'^\u— (ei + ff) — {e2 + ff),v + (ei + ff) + {("2 + v)) . 
Analogously, we can construct a surjective homomorphism from a subgroup of h\^''^\u, v) onto 
Hi^-'^\u-{ei + rf) - (£2 + ??), « + (?i + ^) + (?2 + rf)). This implies that dr iff '^^) < 

El + £2 + 2?7. Since 77 can be taken arbitrari 
so proving the triangle inequality for dr- 



£1 + 62 + 277. Since 77 can be taken arbitrarily small, it follows that dr ( h\^''^\ \ < £1+62, 



□ 

As a simple example, we can easily check that d^ takes the value 1 between the first two persistent 
homology groups represented in Figure [TJ and the value 00 between the first and the third and between 
the second and the third persistent homology groups represented in the same figure. 

Remark 2.3. We observe that the pseudo-distance dr is not a distance. For example, we can consider the 
two Abelian groups Z (obtained by adding infinite copies of the integer numbers) and ^20^^]^ Z. 

We can find two topological spaces X, Y such that their homology groups in degree 1 are these two 
groups, respectively. If we take the two filtering functions ip : X ^ M., ij; : Y ~> M. with ip = ip = 0, 
we have that iff '"^^ (u, v) = Z and iif ''^^ (w, w) = Z2 © Z for any (u, v) G A+ with u > 0. 

Obviously, if (u, v) G A+ with w < we get that both ijf (m, v) and H^^'^^ {u, v) are the trivial 
group. In any case, it is easy to find a surjective homomorphism from ijf (u, v) onto H^'^^ [u, v) 
and a surjective homomorphism from Hj^'^\u,v) onto Hj^^'''^\u,v), so that dr (iif '"^^ iif ''''') = 0. 



However, the groups Hj.^''^\u,v), Hj^'^\u,v) are not isomorphic, for {u,v) e A+ with u > 0. 

2.1. Relationship betvi^een dx and dmatch in the 1-dimensional case, if there is no torsion. 

Let us go back to the case of filtering functions taking values in M. In this case, if the Abelian group 
G of homology coefiicients is also a field, it is well known that the persistent homology groups can be 
described by persistence diagrams, under the assumption that the topological space is triangulable (cf. 
|11| for tame filtering functions, and [5] for continuous filtering functions). The distance that is usually 
used to compare persistence diagrams is the matching distance dmatch- 
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In plain words, the persistence diagram is the collection of all points (u, v) where u and v represent 
the "birth" and "death" of a homology class, varying the sub-level sets of the topological space X with 
respect to the filtering function (p : X ^ R. Each one of these pairs (m, v) is endowed with a multiplicity. 
By definition, all the pairs (u, u) belong to the persistent diagram, each of them counted with infinite 
multiplicity. The matching distance dmatch between two persistence diagrams is the infimum of the 
cost of a matching (i.e. a bijective correspondence) between two persistence diagrams. The cost of 
a matching / is defined to be the maximum displacement of the points of the persistence diagrams, 
induced by /. Each displacement is measured by the max-norm on the real plane. 

For a formal and precise definition of persistence diagram and dmatch, we refer the interested reader 

to [nug. 

It is natural to wonder if there is any relationship between dx and dmatch , under the hypotheses that 
we have just made. 

We can prove the following result. 

Theorem 2.4. Let X, Y be two triangulable spaces endowed with two continuous functions tp : X ^ R, 
: Y ^ M.. Let us consider the persistent homology groups h\^'^\ H^'^^ with homology coefficients 
in a field. Let and D^, be the persistent diagrams associated with the filtering functions tp and ip, 

respectively. Then, for any integer k, dr (^h[^''^\ Ilj^'^^^ < d„iatch{D<4>,D^). 

Proof. Since the coefficients are taken in a field, our persistent homology groups have no torsion, and 
they can be described by their persistent Betti numbers functions j3^, j3^, where we set l3,p{u,v) = 
dimHj^'^''^\u,v), f3^{u,v) = dimHj^'^\u,v), for {u,v) £ A+. Obviously, for each fc G Z, we have 
different persistent Betti numbers functions (5^ of ip (which should be denoted /3<^,fc, say) but, for the 
sake of notational simplicity, we omit adding any reference to k. 

Let us assume that the matching distance dmatch between our persistence diagrams takes the value 
A. Then, from the Representation Theorem 2.17 in [5] (cf. also the k-triangle Lemma in it follows 
that for every A > A and every pair {u,v) G A+ the inequalities I3^{u — X,v + X) < l3^{u,v) and 
/3^(u — A, w + A) < I3ip{u, v) hold. This fact immediately implies that for each {u, v) E A+ 

(1) a surjective homomorphism from a subgroup of h'"^''^^ (u, v) onto h'^'^^ (u — A, u + A)} exists; 

(2) a surjective homomorphism from a subgroup of Hj^''''\u,v) onto iJp'^'(u- A,i; + A)} exists. 

Hence dx {^L[\^''^\h^'^^^ < A for any A > A. This implies our thesis. □ 

Theorem |2.4| leaves open the possibility that dr and dmatch coincide when the persistent homology 
groups have no torsion and hence dmatch is defined. This does not happen. Indeed, for any natural 
number to > 2 we can give an example for which d-r {^L[\^'^\ H^^'^^^ < ^ • dmatch{D^p, D^): 

Example 2.5. Let us consider a natural number to > 2. It is easy to define two triangulable spaces 
endowed with filtering functions whose persistent diagrams consist of the points (0, oo), (0, to), (1, m -|- 
1), . . . , (i, TO+i), . . . , (m, 2to) and (0, oo), (0, to), (1, to+1), . . . , (i, to+z), . . . , (to— 1, 2to— 1), respectively. 
Here all the points have multiplicity one. In Figure [2] we can see the case to = 4. The values of the 
corresponding persistent Betti number functions are displayed in Figure |3] It is easy to check that in 
this case dr and dmatch take the value 1 and y, respectively. 

Theorem |2.4| and Example |2.5| show that the pseudo- distance dr is strictly less informative than the 
matching distance dmatch- However, we recall that dr can be applied also in presence of torsion, while 
dmatch cannot. 

2.2. Stability of the pseudo-distance dx- Another relevant reason to study the pseudo-distance dx 
is its stability. This stability can be expressed as a lower bound for the natural pseudo-distance 5. 
Before proceeding, let us recall the definition of 5. 
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Figure 2. It is easy to check that the matchmg distance dmatch between the two 
persistent diagrams displayed in this figure equals 2. 




Figure 3. Values of the persistent Betti numbers functions represented by the persis- 
tent diagrams in Figure [2j In the picture on the right, the dotted triangle represents 
the part of domain where the persistent Betti numbers functions differ from each other. 
It is easy to check that the pseudo-distance cIt between the persistent homology groups 
described by the two persistent diagrams displayed in Figure [2] equals 1. 



Definition 2.6 (Natural pseudo-distance). Let X^Y be two compact topological spaces endowed with 
two continuous functions 9? : X — > M", -0 : y — > M". The natural pseudo-distance between the pairs 



{X, 0) and (Y, denoted by ,5 ((X, {¥, V^)) , 



is 
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(i) the number inf/i \\0 — "0 ° ^||oo where h varies in the set of all the homeomorphisms between X 
and Y , li X and Y are homeomorphic; 

(ii) oo/ii X and Y are not homeomorphic. 

We point out that the natural pseudo-distance is not a distance because it can vanish on two distinct 
pairs. However, it is symmetric, satisfies the triangular inequality, and vanishes on two equal pairs. 

The natural pseudo-distance has been studied in [131 HH [IZj in the case of scalar-valued filtering 
functions on manifolds, and in |22j in the case of vector-valued filtering functions on manifolds. 

We are now ready to show that the pseudo-distance dr is stable. We begin by proving a useful 
lemma. 

Lemma 2.7. Let X,Y be two homeomorphic compact topological spaces endowed with two continuous 
functions if : X ^ M", ip : Y M". Assume that h : X ^ Y is a homeomorphism such that 
— -0 o h\\oo < e. Then for every pair {u,v) G A+ a subgroup S of H^''^\u,v) and a surjective 
homomorphism / : S" — > H^^'^'^ [u ^ e,v + e) exist. 

Proof. The homeomorphism h induces a chain isomorphism h^. Obviously, takes cycles to cycles 
and boundaries to boundaries. Let us consider the set S of the elements of H^'^\u,v) that contain 
at least one cycle whose image by /i* is a cycle in ^(-0 < u — e). We claim that 5 is a subgroup of 
h\^''^\u,v). Indeed, if ai,a2 G S then two cycles 71 e ai, 72 G 02 exist, such that /i* (71), /i* (72) are 
cycles in y(0 <u — e). Hence the cycle /i*(7i) -f /i*(72) — ^*(7i + 72) is a cycle in Y{tp ^ m — e), too. 
Since 71 -I- 72 G oi -|- 02, it follows that the equivalence class ai -I- 02 belongs to S. That means that S 
is closed with respect to the sum. 

Furthermore, if a G 5 then a cycle 7 G a exists, such that /i*(7) is a cycle in Y{il) < u — e). Hence the 
cycle ft.* (—7) = —ft* (7) is a cycle in Y{'ip ^ u — e), too. Since —7 G —a, it follows that the equivalence 
class —a belongs to 5*. That means that S is closed with respect to the computation of the opposite. 

In summary, S' is a subgroup of h'^'^\u^v). 

Now, let us define the homomorphism f : S ^ h'^'^\u ~ e,v + e). In the following, for each a E S, 
let us fix a cycle 7a G a such that h^,{ja) is a cycle in Y{'ip ^ w — e). For every a £ S we define f{a) as 

the element of Hj^'^'\u — e,v + e) that contains ft* {7a) (i.e., in symbols, [ft*(7a)])- 

We claim that / is a homomorphism. In order to do that, let us consider three elements ai, 02, a G 
h\^''^\u, v), such that ai + 02 — a. Obviously, 

/(fli) +/(a2) = [ft*(7ai)] + [Kha2)] = [h*ilaj + Kila^)] = [h*{la^ +702)] 

/(a) = [ft* (7a)]. 

Now, [ft*(7ai+7a2)]-[^*(7a)] = [^i* (7ai +7a2 ) ^ (7a )] = [K{lat+lao_-la)\- Sinceai-|-a2 = a, a(fc-|-l)- 
cycle T in X{lp ^ v) exists, such that dr = 7^1 +7a2 ~la- Given that — ■0oft||oQ < e, ft*(T) is a (fc+1)- 
cycle va.Y {tp < v-\-€) . Moreover, recalling that ft* is a chain map, ^h^,{T) = h^,{dT) = ft*(7ai +7a2 ^la), 
so that ft*(7ai + 7a2 ^ 7a) is a boundary in :< v + e) and hence [ft*(7ai + 702 ~ 7a)] is the null 
element in H^^''^\u ~ e,v + e). This proves that /(oi) -I- /(a2) — f{a). 

Finally, we claim that / is surjective. Indeed, let b G H^^'^'^ {u — e,v + e) and /3 G 6, so that is a 
cycle in :< u — e). Let a be the cycle ((ft*)""'^) (/3) = (ft~^)* (/3). Given that \\ip o h^^ — ipWoo — 
\\<fi — Tp o h\\ao < e, a is a cycle in X{ip ^ u). Therefore the element a G h\^''^\u, v) the cycle a belongs 
to is an element of the set S. Then /(a) = [ft*(7a)]- Now, 7q and a arc homologous in X{0 < v). Once 
again because of the inequality \\lp — ijj o h\\ao < e, it follows that ft* (7a) and (3 = h^,{a) are homologous 
in Y{ip :< V + e). Hence /(a) — b, and / is proven to be surjective. 

□ 
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Now we can prove that dx gives a lower bound for the natural pseudo-distance S. 

Theorem 2.8. Let X, Y be two compact topological spaces endowed with two continuous functions 
if:X^W\i^:Y^W\ Then dr i?^^''^^^) < 5 ((X, {¥, . 



Proof. Let us set 5 ^ S \^{X,(p), {Y,->p)j. Because of the definition of natural pseudodistance, for each 
> we can find a homeomorphism h^i : X ^ Y such that ||<^— V'°/ii)||oo < S+rj. By applying Lemma 2.7 
to both the homeomorphisms /i,, : X — > y and h^^ : Y X , we get that dx {^H^'^\ nj^'^''^ < S + r]. 
Since rj can be chosen arbitrarily close to 0, our thesis follows. □ 



Previous Theorem |2.8| gives an easy way to get lower bounds for the natural pseudo-distance 6: 

Corollary 2.9. Let X, Y be two compact topological spaces endowed with two continuous functions 
ip : X ^ M", ijj -.Y M". // no surjective homomorphism exists from a subgroup of H^'^\u,v) onto 

Hj^'"^^ {u', tT)}, then S ((X, 0), (Y, iffj > min, min{wi - uj, v'^ - v,}. 

Proof. It follows from Theorem |2.8[ observing that our hypotheses easily imply that the inequality 
dr (h[^'^^ , -ffi,^'''^^) > min,: niin{ui - u'^, v[ - holds. □ 

Our last result illustrates the relationship between dx and the change of the filtering function, mea- 
sured by the max-norm. It shows that d^ is a stable distance. 

Corollary 2.10. Let X be a compact topological space endowed with two continuous functions : X ^ 
M", if:X^W^. Then dr [h';^'^\h[^''''^^ < - i/'lloo- 



Proof. It follows from Theorem 



2.8 



observing that (5 < ||(^- V^'lU- □ 



Conclusions 



The main contribution of this paper is the proof that it is possible to compare persistent homology 
groups in a stable way also in presence of torsion, although no complete representation by persistent 
diagrams or other compact descriptors is available, at the time we are writing. Our approach, based on 
the new pseudo-distance dx, opens the way to experimentation in the case that the homology coefficients 
are taken in an Abelian group instead of a field. In this case, the classical matching distance between 
persistent diagrams cannot be applied without forgetting the information contained in the torsion of 
the persistent homology groups. We plan to extend these ideas to other concepts and problems in 
Multidimensional Persistence. 



Acknowledgments. Research partially supported by DISTEF. The author thanks Sara, and all his 
friends in the "Coniglietti" Band. 



References 

[1] S. Biasotti, A. Cerri, P. Frosini, D. Giorgi and C. Landi, Multidimensional size functions for shape comparison, 
Journal of Mathematical Imaging and Vision, vol. 32, no. 2, 161—179 (2008). 

[2] F. Cagliari, B. Di Fabio, M. Ferri, One-dimensional reduction of multidimensional persistent homology, Proc. Amer. 
Math. Soc, vol. 138, no. 8, 3003-3017 (2010). 

[3] F. Cagliari, C. Landi, Finiteness of rank invariants of multidimensional persistent homology groups. Applied Mathe- 
matics Letters, in press (available online). 

[4] G. Carlsson, S. Gurjeet, A. Zomorodian, Computing Multidimensional Persistence, ISAAC '09: Proceedings of the 
20th International Symposium on Algorithms and Computation, 730-739 (2009). 



10 



p. FROSINI 



G. Carlsson, A. Zomorodian, A. Collins, L. J. Guibas, Persistence Barcodes for Shapes, International Journal of 
Shape Modeling, vol. 11, no. 2, 149-187 (2005). 

G. Carlsson, A. Zomorodian, The theory of multidimensional persistence, SCG '07: Proceedings of the twenty-third 
annual symposium on Computational geometry, Gyeongju, South Korea, 184—193 (2007). 

G. Carlsson, A. Zomorodian, The Theory of Multidimensional Persistence, Discrete and Computational Geometry, 
vol. 42, no. 1, 71-93 (2009). 

A. Cerri, B. Di Fabio, M. Ferri, P. Frosini, and C. Landi. Multidimensional persistent homology is stable. Technical 
Report 2603, Universita di Bologna, 2009. available at http://amsacta.cib.unibo.it/2603/ 

A. Cerri, M. Ferri, D. Giorgi, Retrieval of trademark images by means of size functions. Graph. Models, vol. 68, no. 
5, 451-471 (2006). 

F. Chazal, D. Cohen-Steiner, M. Glisse, L.J. Guibas, S.Y. Oudot, Proximity of persistence modules and their diagrams, 
SCG '09: Proceedings of the 25th annual symposium on Computational geometry, Aarhus, Denmark, 237-246 (2009). 
D. Cohen-Steiner, H. Edelsbrunner, J. Harer, Stability of Persistence Diagrams, Discrete and Computational Geom- 
etry, vol. 37, no. 1, 103-120 (2007). 

D. Cohen-Steiner, H. Edelsbrunner, J. Harer, Y. Mileyko, Lipschitz Functions Have Lp-Stable Persistence, Founda- 
tions of Computational Mathematics, vol. 10, no. 2, 127-139 (2010). 

M. d'Amico, P. Frosini, C. Landi, Natural pseudo- distance and optimal matching between reduced size functions. 
Acta AppUcandae Mathematicae, vol. 109, no. 2, 527-554 (2010). 

V. de Silva, R. Christ, Coverage in sensor networks via persistent homology, Algebr. Geom. Topol., vol. 7, 339-358 
(2007). 

P. Donatini, P. Frosini, Natural pseudodistances between closed manifolds. Forum Mathematicum, vol. 16, no. 5, 
695-715 (2004). 

P. Donatini, P. Frosini, Natural pseudodistances between closed surfaces. Journal of the European Mathematical 
Society, vol. 9, no. 2, 231-253 (2007). 

P. Donatini, P. Frosini, Natural pseudodistances between closed curves. Forum Mathematicum, vol. 21, no. 6, 981-999 
(2009). 

H. Edelsbrunner, J. Harer, Persistent homology — a survey, Contemp. Math., vol. 453, 257-282 (2008). 
H. Edelsbrunner, J. Harer, Computational Topology: An Introduction, American Mathematical Society, 2009. 
H. Edelsbrunner, D. Letscher, A. Zomorodian, Topological persistence and simplification. Discrete & Computational 
Geometry, vol. 28, no. 4, 511-533 (2002). 

M. Farber, Invitation to topological robotics, European Mathematical Society, Zurich Lectures in Advanced Mathe- 
matics, 2008. 

P. Frosini, M. Mulazzani, Size homotopy groups for computation of natural size distances. Bulletin of the Belgian 
Mathematical Society, vol. 6, no. 3, 455-464 (1999). 

R. Christ, Barcodes: the persistent topology of data. Bull. Amer. Math. Soc. (N.S.), vol. 45, no. 1, 61-75 (2008). 
D. Moroni, M. Salvetti, O. Salvetti, Multi-scale Representation and Persistency for Shape Description, MDA '08: 
Proceedings of the 3rd international conference on Advances in Mass Data Analysis of Images and Signals in Medicine, 
Biotechnology, Chemistry and Food Industry, Leipzig, Germany, 123-138 (2008). 

A. Verri, C. Uras, P. Frosini, and M. Ferri, On the use of size functions for shape analysis, Biol. Cybern., vol. 70, 
99-107 (1993). 

A. Zomorodian, Topology for computing, Cambridge Monographs on Applied and Computational Mathematics 16, 
Cambridge University Press, 2005. 

Patrizio Frosini, ARCES, Universita di Bologna, via Toffano 2/2, 1-40135 Bologna, Italia 

DiPARTIMENTO DI MATEMATICA, UnIVERSITA DI BOLOGNA, P.ZZA DI PORTA S. DONATO 5, 1-40126 BOLOGNA, ItALIA, TEL. 

+39-051-2094478, fax. +39-051-2094490 
E-mail address: frosiniOdm.unibo.it 



